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Abstract: In this paper, a fiber optic based sensor capable of fault detection in both radial 
and network overhead transmission power line systems is investigated. Bragg wavelength 
shift is used to measure the fault current and detect fault in power systems. Magnetic fields 
generated by currents in the overhead transmission lines cause a strain in magneto strictive 
material which is then detected by Fiber Bragg Grating (FBG). The Fiber Bragg 
interrogator senses the reflected FBG signals, and the Bragg wavelength shift is calculated 
and the signals are processed. A broadband light source in the control room scans the shift 
in the reflected signal. Any surge in the magnetic field relates to an increased fault current 
at a certain location. Also, fault location can be precisely defined with an artificial neural 
network (ANN) algorithm. This algorithm can be easily coordinated with other protective 
devices. It is shown that the faults in the overhead transmission line cause a detectable 
wavelength shift on the reflected signal of FBG and can be used to detect and classify 
different kind of faults. The proposed method has been extensively tested by simulation 
and results confirm that the proposed scheme is able to detect different kinds of fault in 
both radial and network system. 

Keywords: current measurement; current transformers; optical fiber; magneto strictive 
devices; power system protection 
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1. Introduction 

Current transformers (CTs) and potential transformers (PTs) are widely used in monitoring power 
systems by sending fault current/voltage information to relay and control rooms at substations. If 
current/voltage in a particular line is out of a pre-set range, relays will send a trip signal to breakers to 
trip. If a primary relay fails to operate and clear a fault, based on pre-set rules, a back up breaker will 
clear the fault. Conventional CTs are iron based which subject to saturation and hysteresis. Most relays 
make decisions and send trip signals based on the root mean square (rms) value of fault current 
detected by the CT. Saturation in the CTs cause the rms value of fault current sensed by the CT to be 
much smaller than the actual value and it can prevent relays from tripping and eventually cause 
instability in the system. 

Optical current transformers (OCT) are becoming more popular in power systems. An OCT can 
offer a better transient response, better accuracy, and wider bandwidth in comparison to traditional CTs 
due to the OCTs' lack of iron core [1]. OCTs are light, small, less expensive, and immune to 
electromagnetic interference (EMI) [2] . Already, a number of Faraday effect current sensors have been 
investigated and successfully implemented [3,4]. 

Although OCTs can measure high current without saturation, lack of proper fault detection and 
classification algorithm prevents them to apply in power systems. Previous researchers have discussed 
fault effect in the FBG wavelength shift [5-7]. Authors here, for the first time, use wavelet transforms 
and ANN algorithm to detect and classify the faults based on FBG wavelength shift signals. This paper, 
proposes an FBG based OCT which can replace conventional CTS. A unique ANN algorithm is used 
to first detect the fault based on FBG wavelength shift and then classify it. The proposed OCT does not 
need any CT or PT for biasing and the proposed ANN algorithm works in both radial and network 
systems. 

2. Magnetostrictive Materials 

Magnetostrictive materials are the part of ferromagnetic materials that transform from one shape to 
another in the presence of the magnetic field. Magnetic field causes internal strain in the 
magnetostrictive material with consequences of expansion of the material in the magnetic field 
direction. In magnetostrictive materials, magnetic field strength is proportional to the square of applied 
strain until eventually the magnetic saturation achieved. Since the basis of expansion is molecular, the 
magnetostrictive materials are very sensitive to strain and have a very fast response [8,9]. Also, due to 
the change in the crystal structure of the material, measurement is repeatable with in milliseconds. 
Among these materials, Terfenol-D, Tb 03 Dy n7 Fe ]9S , an alloy of Terbium, Dysprosium and Iron, has the 

highest strain in magnetic field. At room temperature, Terfenol-D can produce about 1,000 ppm which 
is large enough to apply to FBG strain sensor. Previously, Sun and Zheng [10] have shown that the 
highest sensitivity in the Terfenol-D in the magnetic field up to 20 kA/m can be achieved with 6.9 Mpa 
prestress. Due to the nature of all giant magnetostrictive materials, applying prestress can cause a better 
sensitivity. However, their response is roughly proportional to the strength of the magnetic field. 
Terfenol-D can be polarized by using a DC biasing field [8,9]. Performance of the Terfenol-D depends 
on the prestress and the DC bias magnetic field. 
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Since AC magnetic field is measured in this experiment, the DC biasing is necessary to shift the AC 
wavelength and prevent changing the polarity of the output while the input is changing. DC biasing can 
be achieved with serial or parallel permanent magnet, or the DC biasing coil. In this experiment, the 
DC biasing method is used due to its simplicity to change the DC magnetic field. In general, biasing 
point of the magnetostrictive material is defined based on applied prestress and DC biasing point. 
Figure l(a,b) show Terfenol-D with and without DC biasing field respectively. As shown in 
Figure 1(a), the DC biasing field, H Q , is chosen such that the slope of the curve is in its maximum 

point. Behavior of magnetostrictive materials is a nonlinear relation which is already described in 
detail [10,11]. The structural design of FBG sensor using Terfenol-D has been shown in the 
Figure 2(b). In this experiment, DC solenoid with 2,600 amp-turn produces the DC biasing field and 
causes wavelength shift in FBG strain sensor attached to Terfenol-D material. The Hysteresis and eddy 
losses are present in the giant magnetostrictive materials [12-14] and they are considered in the sensor 
model. 

Figure 1. Biasing a Terfenol-D (a) with the DC biasing; (b) without the DC biasing. 



3. FBG Sensors 

FBG sensors can work as arrays for real time measurement of temperature, strain, and pressure in 
the systems. Optical fiber sensors have numerous advantages such as electrically passive operation, 
EMI immunity, high sensitivity, and multiplexing capabilities which make them a perfect candidate to 
use in power systems. FBG sensors are commonly used for strain and temperature measurement and 
they can measure strain up to +5,000 fie and temperature ranges from -40 °C to +120 °C Strain causes 
change in the grating pitch and the fiber index of the sensor. The sensed strain in FBG sensor is then 
coded directly into the wavelength and can be detected as wavelength shift. FBGs reflect a narrowband 
of light and transmit all other wavelengths. In other words, FBG is an optical fiber that works as a filter 
for a particular wavelength. The principal of a FBG based sensor is to detect the reflected Bragg 
wavelength shift due to changes in temperature, strain, or pressure. The Bragg wavelength is defined as 
follows [15]. 




(a) 



(b) 



K = 2n eff A 



(1) 
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where, A is the grating pitch and n eff is the effective index of the fiber core. Bandwidth of the Bragg 

reflected signal depends mainly on grating length and it is typically around 0.05 to 3 nm. 

OCT used in this experiment consists of FBG strain and temperature sensors and a giant 

magnetostrictive material bounded together with epoxy. As described earlier, biasing point of the 

sensor can be adjusted by mechanical prestress on magnetostrictive material and the amp-turn of the 

DC biasing solenoid. One of the main advantages of the OCT sensor, which can substantially reduce 

the cost of implementing power system, is that the OCT sensor does not require any CTs and PTs for 
biasing to detect the fault in the system. Terfenol-D (Tb x Dy x _ x Fe y ) was chosen as the magnetostrictive 

material due to its fast response (<0.1 ms) [16,17] and capability of producing a substantial amount of 
strain in the magnetic field. A broad band source was used to illuminate the FBG through a single 
mode fiber optic and optical coupler. The principle of FBG is well known [15] and will not be 
explained in detail here. A small change in strain and temperature both cause wavelength shift on the 
reflected FBG wavelength by following equations: 

^ = ks + a„ Temp AT. (2) 
A) 

Sn 

a T = JL . (3) 

where, A/tis the wavelength shift, ^is base wavelength in nominal temperature and strain, AT is the 
temperature change, k is the gage factor which is normally 0.78, cc nTemp is the change in index of 

refraction with respect to temperature, and e is strain. In this paper, the following parameters were 
measured through experiment: 2 0 =1528 nm, AA = 1.17 pm /(^strain) ,and wavelength shift is 

AA = 13.48 pm/°C. Also, separate FBG temperature sensor with 27.1 mm probe, ±\.lpml°C sensitivity 

and 0.3 s response time is used to measure temperature and compensate temperate effect. The system is 

calibrated for ambient temperature of 25° degree. 

4. Development of the Concept 

In each phase of the overhead transmission line one sensor is installed and all the sensors should be 
placed as close as possible to its own conductor as shown in the Figure 2(a). They should be placed as 
far as possible from the remaining conductors. Detailed design of sensors and their placement have 
been shown in the Figure 2(b). Since magnetic field is inversely proportional to the distance squared, 
the effect of other phases in this sensor is so small they can be easily ignored in the simulations. AC 
current measurement with FBG has been extensively discussed earlier [18-23]. If overhead 
transmission line is assumed straight infinite conductor, its magnetic field can be expressed as: 

H = — (4) 
2nd 

where, d is the distance from the center of Terfenol-D to the conductor. In our study, two FBG current 
sensors installed at both ends of the transmission line sense strain in both ends of transmission line. 
The transmission line is an assumed bundle with geometer mean diameter (GMD) of 11.07 m and 
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geometric mean ratio (GMR) of 0.127 m with length of 150 km aluminum conductor, steel reinforce 
(ACSR), 60 Hz, and 0.55 per unit impedance. Base values for the system are: 100 MVA and 138 kV. 



Figure 2. (a) Sensor placement in the overhead transmission line, (b) Sensor detail design 
including Terfenol-D, fiber optic cable, and strain sensor. 
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Figure 3. (a) Block diagram of the system; (b) Experimental setup. (DC coil, Terfenol-D 
rod, strain sensor, and temperature sensor). 
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The block diagram of the system and experimental setup are shown in Figure 3(a,b). V reference and 
V fitter are the output of the reference photodiode and filter photodiode, respectively. The reflected light 
(V reference) * s sent through a separate path to the reference photodiode. Previously, Melle and Liu [24] 
have shown that the ratio V fll1er I V reference is proportional to the fiber Bragg wavelength shift: 
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-^ = K{X B -X 0+ ^). (5) 

reference \ 7t 

where, constant K is the slope of the edge or tunable filter, X 0 is the wavelength at which the edge 
filter has zero transmittance, A B is the peak wavelength of the grating, and AA is the spectral width of 
the Bragg reflected spectrum. This method has been chosen to measure the Bragg wavelength shift 
because of its simplicity to install, less complexity and fast response. In this experiment, a DC coil 
made of copper AWG 24 wire with 2,600 amp-turn is used with excitation current of 0.55 A and 
resistance of 18.8 £1. 



5. Simulation Result 

5.1. Radial System 

A typical high voltage power system using electrical transient analyzer program (ETAP) is 
simulated with the assumption of delta-wye connection for all transformers and delta connection for all 
loads in the system. Horizontal configuration is used and the sensors are placed 3 cm away from each 
phase. 

With changing the DC biasing point of the sensor, this method can be applied to different pole 
structures, different cable sizes and cable configurations. A typical radial system with 4 buses is shown 
in Figure 4. In this simulation all equipment is effectively grounded and three sensors are considered 
for each transmission line. 

Figure 4. One line diagram of radial system with sensor location. 



13-8 kV 
Bual 



138 kV 
Bus 2 



138 kV 
Bus3 



6195.4 A 



' 6195.4 



."SI 

-a— 



S3 



• 493 



S4 
— o — 



133 A 
A S2 



S6 



-►372.3 



138 kV 
Bus 4 

S7 

— o- 



a 95 - 8 » • 



6 SS 



8 A 



A SB 



+278.1 A 
A SS 



All the sensors are installed close to the protection relay and are shown with word S followed with 
the sensor number. OCTs are shown by small squares in the system which is installed in all three 
phases and they continuously monitor the magnitude and the phases of the system. 

Increase in the current due to fault in part of the system will cause an increase in the magnetic field 
around the cable. When a short circuit happens in the overhead transmission line, the Bragg 
wavelength of OCT sensed by downstream sensors is shifted due to drastic current change in the 
system. All different kinds of fault between sensor S 3 and ^including single line to ground (SLG), 

double line (DL), and three line to ground have been simulated and the sensors response have been 
shown. The High voltage power system investigated in the paper is assumed to operate at the normal 
condition and the systems are balanced with positive sequence. Since the power system is assumed to 
be balanced and power system analyzing software such as ETAP assume the system in normal 
operation stays balanced, all sensors sense the same wavelength shift for all phases. 
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Figure 5. Fault response at bus 4. (a) System current in the normal condition, (b) Bragg 
wavelength magnitude for each sensor in normal operation of the system, (c) Sensor Bragg 
wavelength magnitude in a three phase to the ground fault at bus 4. (d) Sensor Bragg 
wavelength magnitude in a SLG fault at bus 4. 
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After short circuit calculation for a power system with ETAP, then all signals are sent to Matlab for 
further processing. Normal operation of the system and sensed current are shown in Figure 5(a). 

Multiple sensors can use one fiber optic strand due to the straightforward ability of multiplexing 
several optical current sensors on one optical fiber. Figure 5(b) shows Bragg wavelength shift in 
normal operation of the system. In three phase fault and three phase to ground fault, system current 
stays balanced and it is the same for each phase. As a result, the wavelength shifts detected by all the 
sensors in aforementioned faults stay the same as shown in Figure 5(b). However, in one line to ground 
fault, two line fault, and two line to ground fault current, the system is unbalanced causing different 
wavelength shift for each phase. Figure 5(c) shows when a three phase to ground fault happens in the 
transmission line at bus number 4. Since the system is radial, all upstream feeders and buses should 
sense the high short circuit current. 

The highest surges at sensors S<j and SV hint the closest sensor to fault location are both S<j and S7. 
Figure 5(d) shows the FBG wavelength shift when the single to ground fault (SLG) happens between 
phase A and ground. The highest surge of wavelength shift is in sensor and it therefore shows that 
the fault is the closest to this sensor. Other upstream sensors and their corresponding breakers can work 
as back up in the short circuit fault of the system. 
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Figure 6 shows DL fault between two phases B and C in the same transmission line at bus 4. Like 
all earlier cases, the highest current surge in the system causes the highest wavelength shift resulting at 
the closest sensor to the fault. The highest surge in this system is sensed by and S7 and the current in 
both phases is increased which means DL fault are happened in the system. 

Figure 6. Bragg wavelength shift in a DL fault. 
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5.2. Network System 

In Figure 7 the one line diagram of simple 60 Hz network system with 7 buses has been shown; 
Both 138/13.8 kV transformer in the system considered having unique leakage reactance of 0.06 pu. 
Transmission line is assumed bundle with GMD = 1 1.07 m and GMR = 0.127 m with length of 150 km 
aluminum conductor, steel reinforce (ACSR), 60 Hz with impedance of 0.55 pu. Base values for this 
system are the same as the radial system: 100 MVA and 138 kV. 

Figure 7. One line Diagram of the system with sensors location and their current in normal 
operation of the system. 
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Newton-Raphson load flow method with precision of 0.001 has been used. One line diagram and the 
load flow of the system in normal operation and the Bragg wavelength shift of corresponding sensors 
have been shown in Figure 7, Figure 8 and Figure 9, respectively. In all faults, the magnitudes of Bragg 
wavelength shift detected by the sensors are increased in comparison to either steady state or normal 
operation of the system. Figure 10 indicates three phase fault at bus 4 near sensor S 7 . 

Figure 8. Load flow of system in the normal operation. 
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Figure 9. Bragg wavelength magnitude in the normal operation of the system. 
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Figure 10. Bragg wavelength magnitude in the three phase fault at bus 4. 
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Figure 10 shows the magnitude each sensor exceeds from its normal operation due to increase in the 
current. All current in the approximate of this bus sense huge Bragg wavelength shift, except S3 which 
is connected to the bus feeding motors and in the short circuit motors always act as generators. 
Figure 11 and Figure 12 show magnitude and phase of sensor at single line to ground fault (SLG) at 
bus 4 (phase A fault to the ground at bus 4.) As shown, all sensors connected to phase A of the system 
see a wavelength shift while the other two phases do not sense major wavelength shift. Wavelength 
shifts simply reveal that the concerning faulted phase is phase A of system. 

Figure 11 and Figure 12 reveal that the highest surges happen at sensor 56 and S7 and the second 
highest surge is at sensor 5§and Sn- Since bus 6 is part of the network system, short circuit current can 
flow to the bus 6 from source to the ground through sensors S6, S7, 5§ and Sn- Thus, the fault location 
limits to these sensors and it should be somewhere between of all aforementioned sensors. 

Figure 11. Bragg wavelength magnitude in the SLG fault at bus 4. 
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Figure 12. Sensor phases in the SLG fault at bus 4. 
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Phase reading should be used to detect DL fault in the network system. Figure 13 and Figure 14 
illustrate magnitude and phase of sensor response under DL fault at bus 4, respectively. As expected, 
the highest surges detected by the sensor are in the sensors 56 and S7. The result hints that the fault 
should be at the bus close to these sensors. At DL fault, both phases at the same time can sense the 
same magnitude of surge for both phases A and B. Similar to the previous case, the second highest 
shifts are at the sensors Ss and S12. Based on the topology of the system, the faulted location can be 
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easily recognized and limited to sensors 5g, 57, 5s, and Sn', however, the sensor 5; and S2 sensed the 
higher shift in comparison to the sensor 5§ and S12, they should be ignored since they are not close to 
the reference sensors 5g , S7 nor did their phase reading reveal fault location close to these sensors. All 
sensors around the faulted bus confirm that the current direction is toward the faulted bus. In this 
research high harmonic current sensors mal-function and backup protection in the ANN network were 
not considered which limit the research. Future research should address the harmonics effect at the 
sensor response, sensor mal-function, and sensor noise based on different pole structure and 
configuration. 



Figure 13. Bragg wavelength magnitude in a DL fault at bus 4. 
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Figure 14. Sensor phases in a DL fault at bus 4. 
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5.3. Fault detection and Classification Algorithm 

All fault detection algorithms have two main tasks including detection and classification. 
Previously, numerous algorithms for transmission line fault detection using artificial neural network 
(ANN), wavelet transform (WT) [25], ANN and fuzzy logic [26], ANN and WT [27] have been 
discussed. The fault detection method discussed here, uses FBG wavelength shift data as an input to an 
ANN algorithm. 

In the classification module the voltage and the Bragg wavelength shift samples are recorded. Then 
the samples are normalized by the largest value in the record and the waveforms are resampled for 
1,200 Hz (20 samples/cycle for a fundamental frequency of 60 Hz). A windowing of the samples form 
a set of inputs to the ANN (phase and zero sequence components). A multilayer perceptron (MLP) 



Sensors 2010, 10 



9418 



analyzes the input samples one by one and at the end the final fault classification is defined by the most 
identified fault type. Both voltage and the FBG wavelength shift samples are normalized to their peak 
values at the first cycle in the record, corresponding to the steady state mode of the power system. In 
both high and extra high voltage transmission systems, fault in the system can cause a voltage sag in 
the system and also in some cases some voltage sags induce a large high frequency on both voltage and 
current signals. In addition, this paper has shown that the short circuit in the system can cause 
detectable FBG wavelength shift in the sensor. The fault detection rules are concluded by analysis of 
the FBG wavelength shift in time domain and then in the first decomposition level of discrete wavelet 
transform (DWT) which contains the highest frequency component. Previously, a successful method 
includes a combination of wavelet and ANN for classification, and detection of faults have been 
implemented and tested in the real power network [27]. A similar algorithm is proposed here as shown 
in Figure 15; however, wavelength shift has been used in lieu of current. The fault detection algorithm 
can be summarized as follows: 

(1) If AA pos « AA pre , it is not a fault. 

(2) If £ mx < E* , it is not a fault. 

(3) If £ imx > E" and AA pos < 0. l4AA pre , it is not a fault. 

(4) If ^ E* and \AA pos - AA pre \ < 0. 14 max{ AA pos , M p jit is not a fault. 

(5) If E^ > E* , and none of aforementioned rules have not met, it is a fault. 

Here, the E nm is the maximum energy of the FBG wavelength shift wavelet coefficients, E* is the 

threshold energy defined by the analysis of the wavelet coefficient energy in different system 
operations such as various kinds of fault, and maintenance in part of the system; AX and AX are 

the FBG wavelength shift values on the first and last recording wavelength shift cycles. 

Figure 15. Flowchart of the proposed fault detection method. 
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To calculate the energy of the wavelength shift, a moving window goes through the wavelength shift 
wavelet coefficients shifting one coefficient each time: 

£,=Ik(«)] 2 (6) 

n=\ 

Where, c/, (n)is the nth wavelet coefficient within the z'th window and the window length is n T 
which is defined by n T =n t /2and n s is the number of samples within one cycle of fundamental 

frequency of 60 Hz. For classification purposes like detection module, both voltage and the FBG 
wavelength shift samples are normalized to their peak value at the first cycle in the record. MLP 
network should be trained before implementing the proposed method. To classify all faults correctly, 
the learning data base should contain all varieties and different kinds of fault scenarios to improve 
ANN capability [28,29]. Output of ANN shows which type of fault is related to the input samples and 
in this study binary coding is used for the ANN's outputs. In the case of presence and absence of fault 
the ANN's output will be 1 and 0, respectively. All kinds of faults in binary format have been shown in 
Table 1. 



Table 1. Ann's binary output of the system. 



Fault 


Phase A 


Phase B 


Phase C 


Ground 


type 


(output 1) 


(output 2) 


(output 3) 


(output 4) 


AG 


1 


0 


0 


1 


BG 


0 


1 


0 


1 


CG 


0 


0 


1 


1 


AB 


1 


1 


0 


0 


AC 


1 


0 


1 


0 


BC 


0 


1 


1 


0 


ABG 


1 


1 


0 


1 


ACG 


1 


0 


1 


1 


BCG 


0 


1 


1 


1 


ABC 


1 


1 


1 


0 


No Fault 


0 


0 


0 


0 



Daubechies 4(db4) [30] mother wavelet is used in the simulation and numerous real records 
obtained from different transmission lines have been used to establish the detection rule threshold E* 
and the value is chosen 0.11. To perform the learning process and build the learning data base, faults 
are simulated in the 138 kV line, 150 km long. Different topologies have been investigated and the 
highest and most accurate result has been achieved by using a topology with 40 neurons in the hidden 
layer. Results of the training, validation, and testing have been shown in Table 2. In total, 420 faults in 
bus 4 of radial and then network overhead transmission line were simulated and 99% success were 
achieved in detection and classification of fault in both feeder type as shown in Table 3 and Table 4, 
respectively. 
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Table 2. Simulation data set. 



Variables 


Training 


Validation 


Test 




10-30-50-70-90 






Fault location (km) 


60-110 


80-130 


100-120-140 







Fault type 



AG-BG-CG-AB-BC-AC-ABG-BCG- 
ACG-ABC 



Fault resistance (Q) 



Phase-Phase: 0.5 and 5 
Phase-Ground: 10, 50 and 100 



Table 3. Transmission line fault detection results in the radial system. 



Real system change 


Expected change 
in system 


Number of 
iterations 


Number of 
Success 


Voltage sag 


No fault 


45 


45 


Line de-energization 


No fault 


70 


70 


Normal operation 


No fault 


25 


25 


AG fault 


AG fault 


70 


70 


AB fault 


AB fault 


70 


68 


ABG fault 


ABG fault 


70 


68 


ABC fault 


ABC fault 


70 


70 






420 


416 



Table 4. Transmission line fault detection results in the network system. 



Real system change 


Expected change 
in system 


Number of 
iterations 


Number of 
Success 


Voltage sag 


No fault 


45 


45 


Line de-energization 


No fault 


70 


70 


Normal operation 


No fault 


25 


25 


AG fault 


AG fault 


70 


69 


AB fault 


AB fault 


70 


67 


ABG fault 


ABG fault 


70 


67 


ABC fault 


ABC fault 


70 


70 






420 


413 



4. Conclusions 

A novel protection method in the overhead power transmission line dependent in FBG optical 
current sensor that uses magnetostrictive material is investigated. These FBG current sensors can 
monitor the whole power system simultaneously and they offer better accuracy, and wider bandwidth 
in comparison to conventional CTs. Both radial and network power systems are designed with ETAP, 
and FBG wavelength shifts are processed with MATLAB software. Simulation and signal processing 
results both confirm the effectiveness of the proposed method. Fault detection rules have been defined 
based on ANN algorithm to detect the faulted bus and the rules almost always were successful. 
However, misclassification can occur in the cases that the incident angle, or fault resistance are 
different than what is used during the ANN's learning process. 
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Both the DC biasing point and the prestress should be defined based on the normal current operation 
of the investigated power system. Higher performance can be achieved by precisely defining prestress 
in magneto strictive material, accurately adjusting the DC biasing point of the sensor, and evaluate all 
possible disturbances in the system. Sensors can detect all kinds of fault in both radial and network 
overhead transmission line irrespective of the fault type, fault location, and fault resistance up to 100 
Q.. The main advantages of the FBG current sensors include their ease of installation, immunity to 
EMI, reduced insulation requirement, smaller size, and they are less expensive in comparison to 
conventional CTs. In this research high harmonic current, sensor mal-function, complex feeder 
network and backup protection in the ANN network were not considered which limit the research. 
Future research should address the harmonics effect at the sensor response, sensor mal-function, and 
sensor noise based on different pole structure and configuration. 
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